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Abstract 

We carried out light curve solutions of ten detached eclipsing eccentric binaries observed by Kepler. The 
formal errors of the derived parameters from the light curve solutions are below 1 %. Our results give 
indications that the components of the eccentric binaries (especially those with mass ratios below 0.5) do 
not follow precisely the empirical relations between the stellar parameters derived from the study of circular- 
orbit binaries. We found the following peculiarities of the targets: (a) the components of KIC 9474969 have 
almost the same temperatures while their radii and masses differ by a factor around 2.5; (b) KIC 6949550 
reveals semi-regular light variations with an amplitude of 0.004 and a period around 7 d which are modulated 
by long-term variations; (c) KIC 6220470, KIC 11071207 and KIC 9474969 exhibit tidally induced ’’hump” 
around the periastron. These are the targets with the biggest relative radii of our sample. We derived the 
dependence of the hump amplitude on the relative stellar radii, eccentricity and mass ratio of eccentric 
binary consisting of MS stars. 

Keywords: stars: binaries: eclipsing - stars: fundamental parameters - stars: individual: KIC 6220470, 
KIC 8296467, KIC 6877673, KIC 9658118, KIC 12306808, KIC 5553624, KIC 9474969, KIC 11391181, KIC 
11071207, KIC 6949550 


1 INTRODUCTION 

The most studies on binaries have focused on near¬ 
circular orbits (Pichardo et al. 2005). However, it is 
generally supposed that eccentric binaries are mainly 
produced as a result from the fragmentation (Bonnell & 
Bastien 1992, Bate 1997; Bate & Bonnell 1997). More¬ 
over, the empirical data reveal that the main-sequence 
binary systems typically have eccentric orbits (Duquen- 
noy & Mayor 1991). The probable reason the eccentric 
binaries to be poorly studied is that they are wide stellar 
systems with long periods (Bate et al. 2002) requiring 
prolonged observations. Recently this condition was sat¬ 
isfied by the huge surveys as ROTSE, MACHO, ASAS, 
SuperWASP, etc. The next important step was made 
by the space mission Kepler (Koch et al. 2010). Due to 
its extended and nearly uninterrupted data set above 
thousand detached systems were discovered, consider¬ 
able part of them on eccentric orbits. 

Kepler did not only provide many new eccentric bi¬ 
naries but the unique precision of its observational data 
transformed these binary stars from objects of the ce¬ 
lestial mechanics to an important field of the stellar 
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astrophysics. They became probes for study of the tidal 
phenomena: mechanisms for circularization of the or¬ 
bits and synchronization of the stellar rotation with the 
orbital motion; impermanent mass transfer occurring 
close to the periastron (Sepinsky et al. 2007a, Lajoie & 
Sills 2011); apsidal motion; tidally excited brightening 
and oscillations. 

The theoretical studies reveal that the secular 
changes of the orbital separation and eccentricity could 
be positive or negative (depending on the mass ratio 
and eccentricity), and could occur on timescales ranging 
from a few million years to a few billion years (Sepin¬ 
sky et al. 2007b, 2009). Thus, the widespread assump¬ 
tion for rapid circularization becomes inapplicable, i.e. 
binaries can remain on eccentric orbits for long periods 
of time. Hence, the binary stars on eccentric orbits have 
important evolutional role. 

The eclipsing eccentric binaries (EEBs) with an ap¬ 
sidal motion provide an important observational test of 
the theoretical models of stellar structure and evolution 
(Kopal 1978, Claret & Gimenez 1993). The coefficients 
of internal structure kj are used to describe the external 
potential of a distorted configuration as a function of its 
internal structure in the form of series expansion (Claret 
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& Gimenez 1991). In most cases only k 2 is important 
since higher harmonics in the apsidal motion decrease 
rapidly. In order to compare with observations an av¬ 
erage value k 2 = (cifci 2 -I- C 2 fc 22 )/(ci + C 2 ) of the contri¬ 
bution of star I and 2 (ki 2 and ^ 22 ) is used. The coeffi¬ 
cients Cl and C 2 are known functions of observable stel¬ 
lar and orbital parameters (eccentricity, radii, masses 
and rotational velocities). The individual contributions 
can be separated in two terms, rotational and tidal. In 
most cases tidal terms dominate the apsidal motions of 
binary systems (Claret & Gimenez 1993). Hence, the 
’’apsidal motion test” could be used as a probe of stel¬ 
lar internal structure if precise values of global param¬ 
eters are available (Claret & Gimenez 2010). There are 
important recent studies of this effect based on the ob¬ 
served apsidal motion of double-lined binaries (Lacy et 
al. 2015, Zasche et al. 2014, Bulut et al. 2014, Gar¬ 
cia et al. 2014, Kozyreva & Kusakin 2014, Harmanec 
et al. 2014, Wilson & Van Hamme 2014, Hambleton et 
al. 2013, Wolf et al. 2013, Claret 2012, Zasche 2012, 
Kuznetsov et al. 2011, Claret & Gimenez 2010, Wolf 
et al. 2010, Gimenez & Quintana 1992, Barembaum & 
Etzel 1995, etc.). 

In addition to the classical Newtonian contribution, 
the observed apsidal motion includes term of the Gen¬ 
eral Relativity (Levi-Civita 1937, Gimenez 1985). In¬ 
formation for the apsidal motions of 128 targets in 
our Galaxy can be found in the catalog of Petrova & 
Orlov (1999). The study of EEBs in close galaxies SMC 
and LMC began in the new millenium (Graczyk 2003; 
Michalska & Pigulski 2005; Michalska 2007; Bulut & 
Demircan 2007; North et al. 2010; Zasche & Wolf 2013; 
Zasche et al. 2014). 

Probably the most amazing peculiarities of the bi¬ 
nary stars on eccentric orbits are the tidally excited os¬ 
cillations (harmonics of the orbital period) and bright¬ 
ening around the periastron. They were theoretically 
predicted by Kumar et al. (1995) to explain the ~1 day 
oscillations with amplitude of 0.002 of B-type star or¬ 
biting a neutron star (radio pulsar PSR 0045-7319 with 
e=0.81 and Por6=51 d). Further this analytic theory 
was used for explanation of: (a) the oscillation with fre¬ 
quency of exactly ten times the orbital frequency of the 
slowly pulsating B-star in the binary HD177863 (De 
Cat et al. 2000); (b) the oscillations of the A-type pri¬ 
mary of the eccentric binary HD 209295 (Handler et al. 
2002); (c) the oscillations of the eccentric, short-period 
early-type binary HD 174884 (Maceroni et al. 2009). 

The theory of the tidally excited phenomena further 
was developed by Willems (2003), Zahn (2005), Willems 
& Claret (2005), Willems (2007), Hernandez-Gomez et 
al. (2011), Gundlach & Murphy (2011), Burkart et al. 
(2012), Song et al. (2013), Borkovits et al. (2014), etc. 

Brilliant confirmations of the theoretical predictions 
of Kumar et al. (1995) were discovered by the Kepler 


mission. KOI 54 exhibits two remarkable features: a pe¬ 
riodic brightening spike of 0.7 % occurring at the pe¬ 
riastron and a 0.1 % ’’beat” pattern of pulsations in 
phase with the brightening events (Welsh et al. 2011, 
Burkart et al. 2012). Thompson et al. (2012) discovered 
the next 16 similar objects (most of them noneclips¬ 
ing) in the Kepler archive and called them ’’heartbeat” 
stars due to their shape of light variability reminis¬ 
cent an echocardiogram. Their light curves are differ¬ 
ent: some dim before they brighten, others dim after 
they brighten, and others show distinct W or M shapes. 
Stellar oscillations at harmonics of the heartbeat peri¬ 
ods of some of these targets were also found. 

Answering to the appeal to use the available resources 
of the Kepler database for additional research we un¬ 
dertook study of some types of binary systems from the 
EB catalog (Dimitrov et al. 2012, Kjurkchieva & Dim¬ 
itrov 2015). The goal of this study was to obtain the 
orbits and parameters of ten eccentric binaries based 
on the Kepler data as well as to search for tidally in¬ 
duced phenomena. 


2 SELECTION OF THE TARGETS 

Above two thousands eclipsing binaries (EBs) have been 
identified and included in the Kepler EB catalog (Prsa 
et al. 2011; Slawson et al. 2011), around 1261 of them 
have been initially classified as detached systems. An 
automate fitting of the light curves have been used for 
determination of their ephemerides. 

We reviewed visually the Kepler EB catalog to search 
for detached binaries with eccentric orbits and found 
around 250 targets which phase difference between the 
primary and secondary minimum differed considerably 
from 0.5 and which durations of the two light minima 
were not equal. Most of the found light curves were 
with quite narrow eclipses (due to the long periods). We 
chose to model ten binaries with relatively long eclipses 
(above 0.01 in phase units) allowing precise light curve 
solutions. 

Table 1 presents available information for the tar¬ 
gets (Prsa et al. 2011): orbital period P; Kepler mag¬ 
nitude kmag] mean temperature Tm', width of the pri¬ 
mary eclipse wi (in phase units); width of the secondary 
eclipse W 2 (in phase units); depth of the primary eclipse 
di (in flux units); depth of the secondary eclipse d 2 (in 
flux units). We added to the available target data in 
Table 1 the phases (p 2 of their secondary eclipses (the 
phases (pi of the primary eclipses are 0.0). 

The detailed review of the Kepler data of our targets 
revealed that their light curves did not change during 
the different cycles and observational quarters. This al¬ 
lows to model any data set. We chosen to use several 
consecutive cycles from the middle quarters of each tar¬ 
get. 
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Table 1 Parameters of the targets from the EB catalog (period P, Kepler magnitude kmag, mean temperature Tm, widths Wi and 
depths di of the eclipses) and phases of the secondary eclipses p 2 


Star 

P Id] 

kmag 

Tm 

Wi 

W2 

di 

d2 

<-P2 

KIC 6220470 

8.144136 

13.960 

7377 

0.045 

0.037 

0.091 

0.010 

0.400 

KIC 8296467 

10.33356 

15.177 

5316 

0.024 

0.016 

0.531 

0.326 

0.625 

KIC 6877673 

36.758887 

13.676 

5819 

0.016 

0.013 

0.273 

0.188 

0.439 

KIC 9658118 

24.059995 

14.158 

6225 

0.017 

0.035 

0.473 

0.457 

0.603 

KIC 12306808 

37.878484 

13.265 

5738 

0.011 

0.013 

0.220 

0.169 

0.555 

KIC 5553624 

25.762071 

14.231 

5358 

0.023 

0.007 

0.449 

0.246 

0.420 

KIC 9474969 

21.570506 

12.462 

6085 

0.027 

0.023 

0.133 

0.116 

0.759 

KIC 11391181 

8.617340 

15.257 

5218 

0.017 

0.019 

0.181 

0.106 

0.607 

KIC 11071207 

8.049635 

13.831 

6427 

0.031 

0.040 

0.279 

0.140 

0.621 

KIC 6949550 

7.841067 

15.144 

5719 

0.033 

0.031 

0.355 

0.350 

0.332 


Table 3 Parameters of the best light curve solutions: orbital inclination i, mass ratio q, temperatures Ti, relative radii ri and relative 
luminosities li of the stellar components 


Star 

i 

Q 

Ti [K] 

T2 [K] 

ri 

r2 



KIC 6220470 

89.50 

0.446 

7875 

4382 

0.1014 

0.0286 

0.9918 

0.0082 


±0.01 

±0.001 

±267 

±70 

±0.0010 

±0.0005 



KIC 8296467 

89.58 

0.776 

5726 

5077 

0.0344 

0.0325 

0.6478 

0.5436 


±0.01 

±0.001 

±8 

±6 

±0.0003 

±0.0004 



KIC 6877673 

89.71 

0.798 

6117 

5706 

0.0312 

0.0169 

0.8234 

0.2144 


±0.01 

±0.004 

±29 

±25 

±0.0010 

±0.0010 



KIC 9658118 

89.93 

0.973 

6212 

6193 

0.0431 

0.0396 

0.5445 

0.8365 


±0.01 

±0.002 

±3 

±3 

±0.0002 

±0.0001 



KIC 12306808 

88.668 

0.807 

5772 

5693 

0.0249 

0.0214 

0.5910 

0.6920 


±0.001 

±0.001 

±4 

±3 

±0.0001 

±0.0004 



KIC 5553624 

89.623 

0.682 

6122 

5120 

0.0244 

0.0200 

0.7772 

0.2866 


±0.002 

±0.001 

±8 

±5 

±0.0003 

±0.0005 



KIC 9474969 

87.94 

0.403 

5737 

5574 

0.0679 

0.0258 

0.8876 

0.1266 


±0.01 

±0.002 

±47 

±45 

±0.0001 

±0.0050 



KIC 11391181 

87.42 

0.7728 

5199 

4778 

0.04805 

0.03125 

0.7868 

0.2709 


±0.01 

±0.0014 

±5 

±4 

±0.00030 

±0.00010 



KIC 11071207 

87.47 

0.4304 

6836 

5913 

0.0793 

0.0446 

0.8518 

0.2861 


±0.01 

±0.0020 

±24 

±19 

±0.0007 

±0.0013 



KIC 6949550 

88.45 

0.9859 

5802 

5784 

0.0504 

0.0502 

0.5058 

0.9770 


±0.01 

±0.0018 

±6 

±5 

±0.0008 

±0.0007 





0,88 _ ' ’ __ 

-0,4 -0,2 0,0 0,2 0,4 

Phase 

Figure 1. Top: the folded light curve of KIC 6220470 and its fit; 
Bottom: the corresponding residuals (shifted vertically by differ¬ 
ent number to save space). Color version of this figure is available 
in the online journal. 


3 LIGHT CURVE SOLUTIONS 

We carried out the modeling of the Kepler data by the 
binary modelling package PHOEBE (Prsa & Zwitter 


1.05 



0.45 

0.30 . . . 

-0.2 0.0 0.2 0.4 0.6 

Phase 

Figure 2. Same as Fig. 1 for KIC 8296467 


2005). It is based on the Wilson-Devinney (WD) code 
(Wilson & Devinney 1971, Wilson 1979, Wilson & Van 
Hamme 2004). PHOEBE incorporates all the function¬ 
ality of the WD code but also provides a graphical user 
interface alongside other improvements, including up¬ 
dated filters for the various recent space missions as 
Kepler. That is why PHOEBE is highly appropriate for 
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Figure 3. Same as Fig. 1 for KIC 6877673 Figure 6. Same as Fig. 1 for KIC 5553624 
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Figure 4. Same as Fig. 1 for KIC 9658118 
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Figure 7. Same as Fig. 1 for KIC 9474969 
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Figure 5. Same as Fig. 1 for KIC 12306808 


Figure 8. Same as Fig. 1 for KIC 11391181 
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Figure 9. Same as Fig. 1 for KIC 11071207 



-0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 

Phase 


Figure 10. Same as Fig. 1 for KIC 6949550 


Table 2 The derived orbital parameters of the targets: eccentric¬ 
ity e, periastron angle ui and periastron phase ipper 


Star 

e 

u> (deg) 


KIC 6220470 

0.1879 ± 0.0002 

215.03 ± 0.01 

0.294 

KIC 8296467 

0.2781 ± 0.0001 

314.48 ± 0.01 

0.696 

KIC 6877673 

0.1704 ± 0.0002 

235.26 ± 0.01 

0.369 

KIC 9658118 

0.3723 ± 0.0001 

66.01 ± 0.01 

0.971 

KIC 12306808 

0.1089 ± 0.0001 

36.74 ± 0.01 

0.878 

KIC 5553624 

0.5204 ± 0.0001 

258.12 ± 0.01 

0.411 

KIC 9474969 

0.4165 ± 0.0001 

354.78 ± 0.01 

0.867 

KIC 11391181 

0.1797 ± 0.0001 

18.94 ± 0.01 

0.854 

KIC 11071207 

0.2256 ± 0.0004 

33.00 ± 0.01 

0.897 

KIC 6949550 

0.2661 ± 0.0001 

186.64 ± 0.01 

0.183 


modeling of the precise Kepler data (Hambleton et al. 
2013). 

We used several considerations for the light curve so¬ 
lutions. 

(1) The almost flat out-of-eclipse parts of the ob¬ 
served light curve gave us the grounds to use the de¬ 
tached mode of modeling. 


Figure 11. Sensibility of our light curve solution of KIC 11391181 
(measured by ) to the mass ratio (the rest parameters last fixed 
at their final values) 

(2) The light curve solution of eccentric binaries by 
PHOEBE also by each light curve synthesis software) 
without any guessed values of the eccentricity e and 
periastron angle oj is quite time-consuming task. That 
is why we calculated preliminary values of these orbit 
parameters by the formulae 

TT 

eocoswo =-[(</i2 - <Pi) - 0.5] (1) 

W2 - Wi 

eosma;o =- (2) 

W2 + Wl 

They were obtained as approximations of formulae (9- 
25) and (9-37) of Kopal (1978). We calculated cq and 
ujq by our expressions (1-2) using the values of ip 2 , Wi 
and W 2 from Table 1 ((/ii=0). The obtained values of Cq 
and oiQ were used as input parameters of PHOEBE. 

(3) The mean temperatures Tm of our targets (Table 
1) allowed us to adopt coefficients of gravity brightening 
0.32 and reflection effect 0.5 appropriate for stars with 
convective envelopes. The only exception was the hot 
primary of KIC 6220470 requiring parameters of star 
with radiative envelope (gravity brightening 1.0 and 
reflection effect 1.0). The contribution of the gravity 
brightening effect (and its coefficients) is not consider¬ 
able for detached binaries as our targets. 

(4) Initially the synchronicity parameters were kept 
fixed at values of unity. 

(5) Initially the primary temperature Ti was fixed 
to be equal to the mean target temperature Tm (Table 
1) that has been estimated using dedicated pre-launch 
ground-based optical multi-color photometry plus Two 
Micron All Sky Survey (2MASS) J, K, and H magni¬ 
tudes (Prsa et al. 2011). 

Our procedure of the light curve solutions was quite 
similar to that of Hambleton et al. (2013) and was car¬ 
ried out in several stages. 
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At the very beginning we input some guessed values 
of the secondary temperature T 2 , mass ratio q, orbital 
inclination i and potentials ni _2 and varied only the 
eccentricity e and periastron angle ui around their input 
values Co and ujq to search for the best fit of the phases 
of the eclipses estimated by the value of x^- It turned 
out that the final values e and w differed from the input 
values up to 10 %. This means that the approximated 
formulae (1-2) could be used successfully for calculation 
of the input values of eccentricity and periastron angle 
of eccentric binaries. 

At the second stage we fixed e and w and varied si¬ 
multaneously T 2 , q, i and (and thus relative radii 
^’ 1 , 2 )- We used linear limb-darkening law with limb- 
darkening coefficients corresponding to the stellar tem¬ 
peratures and Kepler photometric system (Claret & 
Bloemen 2011). 

Finally, to adjust the stellar temperatures Ti and T 2 
around the mean value Tm we used the procedure de¬ 
scribed in Dimitrov & Kjurkchieva (2015). 

The final parameters of the eccentric orbits are given 
in Table 2 while Table 3 contains the parameters of the 
stellar configurations. The synthetic curves correspond¬ 
ing to the parameters of our light curve solutions are 
shown in Figs. 1-10 as continuous lines. 

Any surface spots and third lights were not necessary 
to reproduce the photometric data and thus the condi¬ 
tion the synchronicity parameters to be fixed at unity 
turned out without consequence. 

The errors in Tables 2-3 are the formal PHOEBE er¬ 
rors. Most of them are smaller than 1 %, excluding the 
temperatures of the stellar components of KIC 6220470 
whose errors exceed 3 %. We attributed these bigger er¬ 
rors to the considerably shallower eclipses of this target, 
especially the secondary one (see Table 1). The small 
errors of the derived parameters by our light curve so¬ 
lutions are natural consequence of the high precision of 
the Kepler data. 

However, one can see two imperfections of the repro¬ 
ducing of the data. 

The first one are the relative bigger residuals dur¬ 
ing the eclipses (Figs. 1-10). Similar behavior could 
be seen also for other Kepler binaries (Hambleton et 
al. 2013, Lehmann et al. 2013, Maceroni et al. 2014). 
We attributed them to the effects of finite integration 
time studied by Kipping (2010): (i) The large inte¬ 
gration times smear out the light curve signal into a 
broader shape (see figure 1 of Kipping, 2010): the de¬ 
tected ingress and egress durations are bigger than their 
natural values (introducing an additional curvature into 
the eclipse wings) and the apparent positions of the con¬ 
tact points are temporally shifted from their true value; 
(ii) The large integration times smear out the curvature 
of the eclipse bottom caused by the limb-darkening ef¬ 
fect. The two effects are inherent to the long-cadence 
Kepler data of our targets used for modeling. Coughlin 


et al. (2011) also found that the long-cadence (integra¬ 
tion time of 29.43 minutes) Kepler data of systems with 
small sum of relative radii significantly alter the mor¬ 
phological shape of a light curve. 

The second imperfection are the small ripples on the 
flat out-of-eclipse sections of some synthetic light curves 
as well as those on the flat bottom of the total eclipses. 
We attributed them to numerical shortcoming of the 
software for the cases of binaries with small relative 
radii of the components. Such effect (with amplitude 
below 0.1 mmag) was firstly noted by Maceroni et al. 
(2009) and interpreted by the description of stellar sur¬ 
faces with a finite number of elements. The result of 
the surface discretisation seems bigger (with amplitude 
around 0.2 mmag) for our targets which components 
have small relative radii. In fact, this is the reason the 
traditional software for light curve synthesis to be not 
applicable for modeling of exoplanet transits (see for in¬ 
stance PHOEBE scientific reference by A. Prsa, 2011). 

We estimated the global parameters of the tar¬ 
get components (Table 4) using the empirical statisti¬ 
cal relations mass-temperature, radius-temperature and 
luminosity-temperature for MS stars (Boyajian et al. 
2013) which were approximated by combinations of lin¬ 
ear functions. Thus, the absolute parameters in Table 
4 are result of interpolation of empirical relations, not 
derived from observational data, and the listed errors 
are due to the interpolation process. 

4 DISCUSSION 

The analysis of the light curve solutions of the chosen 
ten eccentric binaries led us to several important results. 

(1) We established high sensibility of the solutions 
(measured by x^) to the mass ratio (Fig. 11). This result 
differed from those of Michalska & Pigulski (2004) and 
Terrell & Wilson (2005) who found that for detached 
eclipsing binaries equally good fits can be obtained in a 
very large range of mass ratios. We attributed the sensi¬ 
bility of our solutions to the mass ratio to the following 
reasons: high precision of the Kepler data; very small 
relative radii of the stars of our sample; eccentric orbits 
of our targets. 

(2) The temperatures of the stellar components cor¬ 
respond to spectral type from late A to middle K with 
dominance of G type. This result is expected taking into 
account that Kepler observed mainly solar type stars. 

(3) Our targets have quite big eccentricities (0.11 < 
e < 0.52). This result is due to the selection effect and 
thus does not contradict to the expected distribution of 
e peaked at 0.0-0.1 (Prsa et al. 2008). 

Our targets do not support the expected relation the 
eccentricity to be smaller for the shorter period systems 
(the shorter period orbits undergo stronger tidal forces 
and they circularize more quickly, leaving a dearth of 
highly eccentric short period binaries). Such a trend of 
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Table 4 Masses Mi, radii Ri and luminosities Li of the target components (in solar units) according to the empirical relations. Their 
errors are due to the interpolation process. 


Star 

Ml 

M 2 

Ri 

Ri 


L 2 

KIC 6220470 

1.75±0.09 

0.80±0.01 

1.95±0.08 

0.73±0.09 

10.9±2.2 

0.18±0.02 

KIC 8296467 

l.OOiO.Ol 

0.88±0.01 

1.20±0.01 

1.04±0.01 

1.47±0.01 

0.52±0.01 

KIC 6877673 

1.05±0.01 

0.95±0.01 

1.31±0.01 

1.19±0.01 

2.63±0.09 

1.41±0.04 

KIC 9658118 

1.09±0.01 

1.08±0.01 

1.33±0.01 

1.32±0.01 

3.01±0.01 

3.01±0.01 

KIC 12306808 

0.96±0.01 

0.95±0.01 

1.20±0.01 

1.18±0.01 

1.51±0.01 

1.38±0.01 

KIC 5553624 

1.05±0.01 

0.86±0.01 

1.31±0.01 

1.05±0.01 

2.63±0.02 

0.56±0.01 

KIC 9474969 

1.06±0.02 

0.94±0.01 

1.32±0.01 

1.16±0.02 

1.47±0.09 

1.14±0.06 

KIC 11391181 

0.89±0.01 

0.85±0.01 

1.07±0.01 

0.92±0.01 

0.61±0.01 

0.32±0.01 

KIC 11071207 

1.40±0.01 

0.98±0.01 

1.53±0.01 

1.24±0.01 

5.62±0.11 

1.90±0.04 

KIC 6949550 

0.97±0.01 

0.96±0.01 

1.21±0.01 

1.21±0.01 

1.58±0.01 

1.51±0.01 


increasing of e with the period has been established 
for evolved MS stars in open cluster M67 (Mathieu & 
Mazeh 1988). 

(4) The orbital inclinations of the targets are quite 
near to 90° (Table 3) that is expected for eclipsing sys¬ 
tems with periods above 8 days. But only two targets, 
KIC 6220470 and KIC 9474969, undergo total eclipses. 

(5) We did not find evidences for apsidal motion of 
our targets. The possible reason is the relative short 
duration of the Kepler observations. Typically apsidal 
periods are at least decade long (Michalska & Pigulski 
2005). Moreover, the systems with apsidal motions are 
with the shortest orbital periods or with the largest sum 
of relative radii for a given eccentricity (Michalska 2007) 
but these conditions are not fulfilled for our targets. 

(6) The out-of-eclipse light of the targets is constant 
within 0.2 %. Only KIC 6949550 consisting of two so¬ 
lar type stars reveals semi-regular light variations (Figs. 
12-13) with amplitudes of 0.002-0.004 on timescales of 
order of 7 d (Fig. 12) which are modulated with a pe¬ 
riod around 400 days (Fig. 13). The possibility these 
low-amplitude variations to be detected is due to the un¬ 
precedent precision of the Kepler data. Their frequency 
analysis is object of future study. 

(7) The review of the light curves of our targets from 
different quarters did not exhibit any long-term vari¬ 
ability (excluding KIC 6949550). 

(8) The mass ratios of the targets are within the range 
0.4-1.0. North et al. (2010) obtain that q for detached 
systems is within 0.8-1.1 and for semidetached and con¬ 
tact binaries within 0.4-0.7. Lucy (2006) found an ex¬ 
cess of binary systems with q > 0.95 in our Galaxy. Our 
result implies that probably the distribution of the mass 
ratio for eccentric binaries differs from that for circular- 
orbit systems. 

(9) We established some linear trend between the ra¬ 
tio r 2 /ri and mass ratio q (Fig. 14, left) for the ten 
members of our small sample. 

(10) The ratios T^/Ti of our targets fall in the range 
0.8-1.0 with one exception, KIC 6220470, which value 
is 0.56 (Fig. 14, middle). Opposite discrepancy has been 
found for HD 174884 (Maceroni et al. 2009) which com¬ 


ponents are with equal temperatures irrespective of the 
negligible secondary eclipse. 

(11) The values pst of the mass ratios, obtained by the 
empirical relation temperature-mass of MS stars (Ta¬ 
ble 4), differed from those determined by our PHOEBE 
models q = qpH (Fig. 14, right). The most glaring case 
is KIC 9474969 with values qst=0.88 and qpH=0-4:. 

Figure 14 implies that the components of the eccen¬ 
tric binaries, especially those with q < 0.5, do not follow 
the empirical relations between the global stellar param¬ 
eters derived by study of circular-orbit binaries (and 
with probable dominance of mass ratios around unity). 
In fact only targets KIC 9658118 and KIC 6949550 with 
mass ratio near unity turned out with almost equal val¬ 
ues of the ratios of luminosities and radii of their com¬ 
ponents obtained by the light curve solutions {h/h and 
r^jrx from Table 2) and by the empirical relations for 
MS stars (L 2 /L 1 and R 2 /R 1 from Table 4). The biggest 
discrepancies between the values of these ratios belong 
to KIC 9474969 with the smallest mass ratio. This re¬ 
sult might be considered as some empirical support 
of the conclusion of Sepinsky et al. (2007a) about the 
difference of the Roche geometry of circular-orbit and 
eccentric-orbit binaries. Another possible reason could 
be that the empirical relations for MS stars have been 
derived on the base of binaries with mass ratios around 
unity. 


5 TIDALLY INDUCED HUMPS 

The tidal forces change the stellar shape (tidal bulges) 
and cause brightness variability due to projection of the 
distorted stellar surfaces on the visible plane (Brown et 
al. 2011; Welsh et al. 2011; Morris 1985). It has double¬ 
wave shape (ellipsoidal variations) in the case of circular 
orbits and light increasing around the periastron in the 
case of eccentric orbits. 

Kumar et al. (1995) created an analytic model of 
tidal phenomena in eccentric binary consisting of point 
source (neutron star) and MS star. The shape of the 
corresponding light curve depends on inclination, angle 
of periastron, and eccentricity while its amplitude (of 
order of mmag) depends on the masses of the objects, 
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Figure 12. Short-term semi-regular light variations of KIC 6949550 
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Figure 13. Long-term modulation of the light variations of KIC 6949550 




Figure 14. Dependence of on mass ratio q (left panel), dependence of T 2 /T 1 on mass ratio q (middle panel) and diagram qst - 

qpH (right panel) 


their internal stellar structure and the orbital separa¬ 
tion at the peristron (Kumar et al. 1995). According to 
this model the shape of the light increase is one-peaked 
for i < 30° but becomes two-peaked with central dip 
(which depth and width increase with i) for the bigger 
orbital inclination. 

Thompson et al. (2012) calculated a grid of solu¬ 
tions to the model of Kumar et al. (1995) and found 
that the larger inclinations cause the light curve firstly 
to increase in brightness and then to decrease (KIC 
3547874), or vice versa (KIC 9790355), depending on 
the angle of periastron, while the bigger eccentricity led 
to shorter duration of the heartbeat event. 


Three targets from our sample, KIC 6220470, KIC 
9474969 and KIC 11071207 (Fig. 15), reveal light fea¬ 
tures around the periastron phase (pper (Table 2). The 
analysis of these observed features led us to the follow¬ 
ing results. 

(a) The observed features around the periastron 
phase (Fig. 15) seem as a ’’hump” (brightening). Their 
shapes differ from the expected ones for the big orbital 
inclinations (with central dips, see figs. 2-3 of Kumar 
et al. 1995 and fig. 5 of Thompson et al. 2012). 

(b) The hump duration is shortest for KIC 9474969 
which eccentricity is above 2 times bigger than those of 
the other two ’’humped” targets (see Table 2 and Fig. 
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Figure 15. Tidally induced brightening of KIC 6220470 (left panel), KIC 9474969 (middle panel) and KIC 11071207 (right panel) 


15). This result supports the conclusion of Thompson 
et al. (2012). 

(c) The tops of the humps do not coincide perfectly 
with the periastron phases: those of KIC 6220470 and 
KIC 11071207 slightly precede periastron phase while 
that of KIC 9474969 slightly delays (Fig. 15). These 
small deviations (phase shifts below 0.04) probably are 
due to the irradiation effect which phase contribution 
depends on the periastron angle. 

(d) Humps were found only for those targets from our 
sample for which ri + r 2 > 0.09 (Table 2). This implied 
that the hump amplitude depends strongly on the rela¬ 
tive stellar radii. An additional argument was that the 
hump amplitude was biggest for KIC 6220470 (Fig. 15) 
which value ri -|- r 2 = 0.13 was the biggest one (Table 
3) among the three ’’humped” targets of our sample. 

To derive the dependence of the hump amplitude on 
the relative radii we applied the formalism of Kumar et 
al. (1995) to eccentric binary consisting of two MS stars 
and obtained the expression 


SF 1 r| r? 

— = 4- -I- 4o-. 

F q{l-e)3 (l-e)3 


( 3 ) 


where F is the total flux of the target. In fact, this 
formula gives the amplitude of the ’’clean” tidally in¬ 
duced hump. The irradiation effect, which contribution 
depends on the periastron angle, superposes the tidally 
excited brightening and changes the hump amplitude. 
Moreover, possible tidally-induced pulsations (Willems 
& Aerts 2002) and Doppler boosting (Bloemen et al. 
2011) could cause additional complications. 

The hump amplitudes of our targets calculated by 
the expression (1) were around 2 times bigger than the 
observed ones. Kumar et al. (1995) obtained just the 
same discrepancy for PSR 0045-7319. This led us to the 
supposition about wrong coefficient of their expressions 
(2 instead of 4). 

We should point out that KIC 6949550 has ri -|- 
r2=0.I but does not reveal a hump. We attributed this 


exception from the rule (existence of detectable hump 
for ri -|- r 2 > 0.09) by its semi-regular variations with 
bigger amplitude which blur the possible hump. 

We assume that the obtained limit of 0.09 of ri -|- r 2 
depends on the data precision and can be reduced for 
future space missions. 

(e) All targets with periastron brightening have q ~ 
0.45. We attributed this fact as a consequence of the 
accidental coincidence of the binaries with q cs 0.45 and 
the systems with biggest relative radii because the de¬ 
pendence (1) of SF/F on the relative radii is stronger 
than that on the mass ratio. 

(f) Humps would not be discern for systems which 
periastron angle uj is close to 90° or 270° because the 
periaston phase for these cases will coincide with some 
of the eclipses. Our targets are not such cases (Table 
2) and thus, we have not missed some humps on this 
reason. 

(g) We did not find oscillatory modes of any tar¬ 
get. This result confirms the conclusion of Kumar et 
al. (1995) that the MS stars have small amplitudes of 
oscillation. 


6 CONCLUSION 

This paper presents the results of determination of the 
orbits and fundamental parameters of ten eclipsing bi¬ 
naries with eccentric orbits from the Kepler archive. 
The obtained results allowed us to derive several con¬ 
clusions. 

(1) The formal errors of the derived parameters from 
the light curve solutions are below 1 % (excluding those 
of KIC 6220470 which exceed 3 %). 

(2) Our light curve solutions imply that the compo¬ 
nents of the eccentric binaries (especially those with 
mass ratios below 0.5) do not follow precisely the em¬ 
pirical relations between the stellar parameters derived 
from the study of circular-orbit binaries. 
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(3) KIC 6949550 reveals semi-regular light variations 
with an amplitude around 0.004 and a period around 7 
d which are modulated by long-term variations. 

(4) We found tidally induced light ’’hump” around 
the periastron phase of three targets: KIC 6220470, KIC 
11071207 and KIC 9474969. 

(5) We derived formula describing the amplitude of 
the tidally induced hump of eccentric binary consisting 
of MS stars. It exhibits that the amplitude of these fea¬ 
tures increases strongly with the relative stellar radii. 

(6) Although we did not find evidences of apsidal mo¬ 
tion of our targets in the framework of the relative short 
duration of the Kepler observations our EEBs present 
appropriate targets for future study of this effect due to 
their big eccentricities. 

Obviously, the numerous and exclusive precise Ke¬ 
pler data deserve precise light curve solutions to enrich 
statistics of the binaries with estimated parameters and 
thus to improve the empirical relations between them. 
Moreover, the tidally induced phenomena provide criti¬ 
cal data to constrain the theories of tidal forces in stellar 
binaries (Euller & Lai 2011, Burkart et al. 2012). 
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